Half-metallic zigzag carbon nanotube dots 
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A comprehensive first-principles theoretical study of the electronic properties and half-metallic 
nature of finite zigzag carbon nanotubes is presented. Unlike previous reports, we find that all 
nanotubes studied present a spin-polarized ground state, where opposite spins are localized at the 
two zigzag edges in a long-range antifferomagnetic configuration. Relative stability analysis of the 
different spin states indicate that, for the shorter segments, spin-ordering should be present even 
at room temperature. The energy gap between the highest occupied and the lowest unoccupied 
molecular orbitals of the finite systems is found to be inversely proportional to the nanotubes 
segments length, suggesting a route to control their electronic properties. Similar to the case of 
zigzag graphene nanoribbons, half-metallic behavior is obtained under the influence of an external 
axial electric field. 
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I. INTRODUCTION 

Magnetism in carbon based materials has attracted 
considerable scientific interest in recent years both from 
an experimental 1-17 and a theoretical 18-53 viewpoint. 
While the origin of magnetic ordering in such systems 
is yet to be fully understood, it has been suggested 
that spin polarization may arise from local structural 
defects, 33,36,41,53 sterically protected carbon radicals 34 
and chemical impurities. 30,40, A unique mechanism for 
spin ordering in graphene based systems is related to 
the appearance of edge states. 18-29,42,45-52 When cutting 
a graphene sheet along its zigzag axis to form a nar- 
row and elongated graphene nanoribbon (GNR), distinct 
electronic states appear, which are localized around the 
exposed edges. 54-59 These states are predicted to carry 
spin polarization, resulting in a well defined magnetic 
ordering. 19-26,28,29,42,45-52 Due to the bipartite hexago- 
nal structure of graphene, the electronic ground state of 
such zigzag graphene nanoribbons (ZZ-GNRs) is char- 
acterized by an antiferromagnetic (AFM) spin ordering, 
where the edge states located at the two zigzag edges of 
the ribbon have opposite spins. Under the influence of 
an external electric field, these systems have been pre- 
dicted to become half-metallic with one spin channel be- 
ing semiconducting and the other metallic, 45,47 thus act- 
ing as perfect spin filters with important implications to 
the field of spintronics. Furthermore, this behavior has 
been shown to be preserved for finite graphene nanorib- 
bons, where the zigzag edge may be as short as a few 
repeating units. 48-52 

Similar to graphene nanoribbons, is has recently been 
suggested that spin ordering may occur on the zigzag 
edges of hydrogen terminated finite sized carbon nan- 
otubes (CNTs). 32,37 Based on density functional theory 
(DFT) calculations within the local spin density approx- 
imation (LSD A), it was predicted that the ground elec- 
tronic state of these systems strongly depends on their 
circumferential dimension ranging from AFM to ferro- 
magnetic (FM) ordering, where both edges of the nan- 
otube bear spins of the same flavor. This conclusion is 
in striking difference from the case of ZZ-GNRs, where 



the ground electronic state was found to have an AFM 
spin ordering regardless of the dimensions of the rib- 

k on 19-26,28,29,42,45-52 

Interestingly, contradicting results have been reported 
in the literature, 33 indicating that the ground state of 
unpassivated zigzag CNTs has a low-spin AFM order- 
ing even for nanotubes that have earlier been considered 
to present a high-spin ferromagnetically ordered ground 
state. 32 

In light of the considerable progress that has been 
made in the synthesis and fabrication of ultra-short 60-64 
and open ended 65 CNTs and ultra-narrow GNRs, 66-69 it 
is desirable to obtain a full understanding of their elec- 
tronic properties. Such understanding may prove impor- 
tant in future applications of these molecules in nano- 
electronic and nanospintronic devices. 



(7,0,4) 

j J J J * ' 

(8,0,4) 

J J J -I 

(9,0,4) 
(10,0,4) 



(7,0,7) 

** »-» *-* * * 

4 i> 4 4 i> 4 
(8,0,7) 



(9,0,7) 

v<KX 

(10,0,7) 



, ****** *4 4~ 

» » • * * <i #•* * * 
(7,0,10) 

(8,0,10) 

4 4 4 4 »« • « 4 4 

** »«; *-4j n »* 

(9,0,10) 

*4 4 4 j -> i j .j 4' 
j j j j j j j j j j 

'4 J -I * * 4 * > * * ' 

(10,0,10) 



(7,0,13) 
J *i 5' i J 5 J ** ** 5* *\ *' ** ' 

j J 4-4 4-4 4-4 * 4 J J J - 

*4_ <M »A **L 

(8,0,13) 

4 4 *-4\ 4-4 4 4 » * 4 4 4~J 

(9,0,13) 

. jj Z ' xj 2 • *J ii xl J . xl 4 . m4 J . xj j 

' 4 j j j j j j j j j j j j 

^* xt-xf »-» 9-» »-» 

(10,0,13) 



FIG. I: The four sets of CNTs studied. The notation 
(N, M, L) indicates a finite segment of a (N, M) CNT with L 
zigzag rings along its principal axis. 

It is the purpose of the present paper to present a com- 
prehensive investigation of spin polarization in hydrogen 
passivated finite zigzag CNTs. Unlike previous reports on 
these systems, 32,37 we find that all zigzag CNTs studied 
present an antiferromagnetically ordered ground state, 
regardless of their diameter and length. Furthermore, as 
is the case with infinitely long 45-47 and finite 52 GNRs, 
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the application of an external axial electric field results 
in half-metallic behavior. 



II. METHOD 

We have studied a set of sixteen finite segments of the 
(7,0), (8,0), (9,0), and (10,0) zigzag CNTs (see Fig. 1). 
For each CNT four segments of 0.93, 1.56, 2.20, and 2.84 
nm in length are considered. All nanotubes considered 
are hydrogen terminated, i.e., each carbon edge atom is 
passivated with a single hydrogen atom. We label the 
CNTs by (N,M,L), where (N, M) is the infinite CNT 
from which the relevant finite segment is derived and L 
is the number of zigzag carbon rings stacked together to 
form the finite CNT (see Fig. 1). This set of nanotubes 
includes all those previously predicted to present a ferro- 
magnetically ordered ground state. 32 

All the calculations presented in this work were car- 
ried out using the development version of the Gaussian 
suite of programs. 70 Spin polarized ground state calcula- 
tions were performed using the screened exchange hy- 
brid density functional of Heyd, Scuseria and Ernzer- 
hof (HSE06), 71-74 which has been tested on a wide va- 
riety of materials and has been shown to accurately re- 
produce experimental bandgaps 75-77 and first and sec- 
ond optical excitation energies in metallic and semicon- 
ducting single walled CNTs. 78 ' 79 The inclusion of short- 
range exact-exchange in the HSE06 functional makes 
it suitable to treat electronic localization effects, 80-84 
which are known to be important in this type of materi- 
als _ is-29,42,45-48,52,56-59 Thig ^ further sup p rted by the 

good agreement, which was recently obtained, between 
predicted bandgaps 85 of narrow graphene nanoribbons 
and measured values. 66-68 To obtain a reliable ordering of 
the different magnetization states, we find it important to 
relax the geometry of the finite CNTs for each spin polar- 
ization. Therefore, unless otherwise stated, all reported 
electronic properties are given for fully optimized struc- 
tures using the double-zeta polarized 6-31G** Gaussian 
basis set. 86 It should be noted that since our calculations 
are performed within a single determinantal framework, 
we can determine only the total spin vector projection 
along a given axis, m s , and not the total spin. This is 
standard in unrestricted Kohn-Sham theory. 



III. RESULTS AND DISCUSSION 

We start by performing ground state calculations on 
the full set of sixteen finite CNTs studied. In Fig. 2 
spin-density maps of the ground state of four representa- 
tive CNT segments are presented. In contradiction with 
previously reported LSDA results, 32,37 which predict the 
(8, 0) CNT derivatives to have an AFM ground state and 
the (7,0) and (10,0) nanotube segments to have a fer- 
romagnetically ordered ground state, our HSE06 calcu- 
lations indicate that all CNTs studied have an antifer- 



romagnetically ordered ground state, with a total spin 
vector projection of m s = and a spatially resolved 
spin density. Similar to the case of finite GNRs, 48,52 one 
zigzag edge of the CNTs segments has a high density of 
spin a electrons (blue color in Fig. 2) while the other edge 
is rich with spin [3 electrons (red color in the figure). 
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FIG. 2: Antiferromagnetic ground state spin density maps 
of the (7,0,7) (upper left panel), (8,0,7) (upper right panel), 
(9,0,7) (lower left panel)), (10,0,7) (lower right panel), finite 
zigzag CNTs as obtained using the HSE06 functional with the 
6-31G** basis set. Red and blue isosurfaces indicate the two 
spin flavors with an isovalue of 0.0015 a^ 3 . 



To understand the origin of this discrepancy between 
the predictions obtained using the semi-local and the 
screened-hybrid functionals, we have repeated the calcu- 
lations using the LSDA functional and the same atomic 
basis set. Interestingly, we encountered considerable con- 
vergence problems when performing the calculations at 
this level of theory. For many of the systems we have 
found it necessary to start the SCF with a pre-designed 
initial guess in order to achieve convergence to the ap- 
propriate spin state. Once convergence was achieved we 
have found that for the (7, 0, 7) and (7, 0, 10) systems 
the LSDA functional does predict a ferromagnetically or- 
dered ground state. Nevertheless, this state is only 5 
meV below the antiferromagnetically ordered state. For 
the rest of the systems considered, whenever convergence 
was achieved, an AFM ground state has been found. 

We attribute this behavior to the self-interaction error 
appearing in the LSDA functional, which tends to con- 
siderably over estimate electron derealization. 87,88 Since 
some of the most important physical characteristics of the 
systems under consideration relate to pronounced local- 
ized edge states, care should be taken when using the 
LSDA functional to capture the correct energetic order- 
ing of the different spin polarized states. As mentioned 
above, the inclusion of Hartree-Fock exchange within hy- 
brid functionals considerably reduces the derealization 
error and thus makes them more appropriate to treat 
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scenarios where electron localization plays an important 
rule in determining the electronic structure of the sys- 
tern. 80 - 84 
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FIG. 3: Energy differences between the antiferromagnetic 
ground state and the above lying higher spin multiplicity fer- 
romagnetic state for the four sets of CNTs segments stud- 
ied, as calculated by the HSE06 functional. feT at room 
temperature is indicated by the dashed magneta line. In- 
set: Energy difference between the antiferromagnetic ground 
state and the non-magnetic closed shell electronic state of the 
(10,0,L) CNT segments. 

To quantify our findings, we study the stability of the 
antifcrromagnetically ordered ground state with respect 
to the above lying ferromagnetically ordered spin state. 
In Fig. 3 the energy differences between the antiferromag- 
netic ground state and the first higher spin multiplicity 
state are presented for the four sets of nanotube segments 
studied. We indicate fc^T (ks being Boltzmann's con- 
stant) at room temperature (T=298 K) by the dashed 
magenta (color online) line in the figure. For the (7, 0, L) 
and (9, 0, L) systems we find the m a — 2 state to be the 
above lying ferromagnetic state. For the (10,0, L) seg- 
ments the first higher spin multiplicity state has m s = 3. 
In the case of the (8, 0, L) segments the higher spin mul- 
tiplicity state changes from m s = 1 for L = 4 to m s = 3 
for L = 7, 10 and 13. As can be seen, the AFM ground 
state of the shorter CNT segments studied is considerably 
more stable than the above lying higher spin state, and 
is expected to be detectable at room temperature. Like 
in the case of GNRs, 45 the energy difference between the 
ground state and the above lying higher spin-state de- 
creases monotonically with increasing CNT length, up to 
a point where the differences become lower than room 
temperature. Interestingly, we find that for the (9, 0) 
segments, which in the limit of infinite length become 
metallic, the decay rate of the energy difference between 
the AFM ground state and the above lying FM state is 
slower than that of the other semiconducting tubes stud- 
ied. Even for this system an exponential fitting of the 



decay curve suggests that at a length of 8 nm the en- 
ergy difference is below room temperature. The inset 
of Fig. 3 shows the energy difference between the AFM 
ground state and the non-magnetic closed-shell electronic 
state of the (10,0, L) segments. The non-magnetic state 
is found to be considerably higher in energy when com- 
pared to higher spin multiplicity states. This, as well, 
is in agreement with the case of GNRs 45,52 and suggests 
that even at room temperature short segments of CNTs 
are expected to present magnetic properties. 
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FIG. 4: HOMO-LUMO gap values as a function of CNT seg- 
ment length for the four sets of CNTs studied as calculated 
by the HSE06 functional. 

Before discussing the effect of an electric field on the 
electronic properties of finite CNTs, it is essential to 
study their ground state characteristics in the absence 
of external perturbations. The length dependence of 
the HOMO-LUMO (highest occupied molecular orbital 
and lowest unoccupied molecular orbital, respectively) 
gap would be the most important parameter to address. 
In Fig. 4 the energy gap as a function of the length of 
the CNTs segments are presented, for the four subsets 
of nanotubes considered. All studied finite zigzag CNT 
segments have sizable HOMO-LUMO gaps, including the 
(9, 0, L) series, which in the limit of L — > oo is expected 
to become metallic. The HOMO-LUMO gap is inversely 
proportional to the length of the CNTs segments for all 
nanotubes studied and its value strongly depends on the 
length, changing by more than 0.8 eV for lengths in the 
range of 0.9 and 2.8 nm. This suggest that careful tai- 
loring of the nanotube length can be used as a sensitive 
control parameter over its electronic properties. 89 

We now turn to check whether the similarity between 
finite zigzag CNTs and finite GNRs extends also to the 
case of half-metallic behavior under the influence of an 
external electric field. 45,52 In Fig. 5 we present the spin- 
polarized HOMO-LUMO gap dependence on an external 
electric field applied parallel to the main axis of the nan- 
otubes (perpendicular to the zigzag edges) for four repre- 
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FIG. 5: Spin-polarized HOMO-LUMO gap dependence on the 
strength of an external axial electric field for the (7, 0, 7) (up- 
per left panel), (8, 0, 7) (upper right panel), (9, 0, 7) (lower left 
panel), and (10, 0, 7) (lower right panel) finite zigzag CNTs as 
calculated by the HSE06 functional. Fixed geometries of the 
relaxed structures in the absence of the external field were 
used. Full and open marks stand for the a and f3 spin gaps, 
respectively. 

sentative finite zigzag CNTs. In agreement with previous 
calculations, 45,47 ' 52 the a and (3 gaps are degenerate in 
the absence of an external field. Upon application of the 
field, electrons having one spin experience an increase in 
the HOMO-LUMO gap while the opposite spin flavor ex- 
periences a decrease in the gap. This gap splitting contin- 
ues up to a point, where one spin channel presents a van- 
ishing gap, thus creating a half-metallic state or a perfect 
spin filter. At this point, due to the increased mobility 
of the metallic electrons, further increase in the exter- 
nal field results in spin transfer between both edges thus 
reducing the total spin polarization and the energy gap 
splitting (clearly seen on the upper left panel of Fig. 5). 
All four representative structures present the same fea- 
tures described above. The main differences observed are 



in the zero field HOMO-LUMO gap and the onset field for 
the appearance of the half-metallic state. Nevertheless, 
they all predict a half-metallic state at an appropriate 
electric field strength. 



IV. CONCLUSIONS 

To summarize, intrigued by reports of contradicting 
results regarding the ground state properties of finite 
zigzag carbon nanotubes, we performed a comprehen- 
sive study of the electronic character of a representa- 
tive set of sixteen CNT segments. Unlike previously re- 
ported results, we found that all finite zigzag CNTs pos- 
sess a spin-polarized ground state with antiferromagnetic 
spin ordering. This state is characterized by high spin- 
density of opposite spins located at the two zigzag edges 
of the molecule. Our HSE06 results predict the antiferro- 
magnetic ordering to be considerably stable with respect 
to higher spin multiplicity states and the non-magnetic 
closed shell state, suggesting that their spin polarization 
is detectable at room temperature. The HOMO-LLIMO 
gap was found to be inversely proportional to the length 
of the zigzag CNT segment. The high sensitivity of the 
gap to changes in the length of the CNT segment sug- 
gests a way to control the electronic properties of such 
systems. Similar to the case of zigzag GNRs, the half- 
metallic nature of finite zigzag CNTs under an external 
in-plane electric field was verified. Due to the recent suc- 
cess of the HSE06 functional in predicting the electronic 
properties of GNRs and CNTs of different nature and 
dimensions, we are confident about the reliability of the 
predictions presented here. 
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